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Abstract: 

Building on earlier work, the dipole subtraction formalism for photonic corrections is 
extended to various photon-fermion splittings where the resulting collinear singularities 
lead to corrections that are enhanced by logarithms of small fermion masses. The differ- 
ence to the earlier treatment of photon radiation is that now no cancellation of final-state 
singularities is assumed, i.e. we allow for non-collinear-safe final-state radiation. More- 
over, we consider collinear fermion production from incoming photons, forward-scattering 
of incoming fermions, and collinearly produced fermion-antifermion pairs. For all cases we 
also provide the corresponding formulas for the phase-space slicing approach, and particle 
polarization is supported for all relevant situations. A comparison of numerical results 
obtained with the proposed subtraction procedure and the slicing method is explicitly 
performed for the sample process e~7 — > e~fi~fi + . 
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1 Introduction 



Present and future collider experiments require precise predictions for particle reac- 
tions, i.e. for most of the relevant processes radiative corrections have to be calculated. 
This task becomes arbitrarily complicated if either the order in perturbation theory (loop 
level) or the number of external particles is increased, or both. In recent years the needed 
techniques and concepts have received an enormous boost from various directions; for a 
brief overview we refer to some recent review articles [1,2]. 

In this paper we focus on real emission corrections involving photons at next-to-leading 
order (NLO). Apart from the integration over a many-particle phase space, here the main 
complication is the proper isolation of the singular parts which originate from soft or 
collinear regions in phase space. To solve this problem at NLO, two different types of 
methods have been developed in the past: phase-space slicing (see, e.g., Ref. [3]) and 
subtraction [4-8] techniques. In the slicing approach the singular regions are cut off 
from phase space in the numerical integration and treated separately. Employing general 
factorization properties of squared amplitudes in the soft or collinear regions, the singular 
integrations can be carried out analytically. In the limit of small cutoff parameters the 
sum of the two contributions reproduces the full phase-space integral. There is a trade-off 
between residual cut dependences and numerical integration errors which increase with 
decreasing slicing cuts; in practice, one is forced to search for a plateau in the integrated 
result within some errors by varying the slicing cut parameters. 

This cumbersome procedure is not necessary within subtraction formalisms which are 
based on the idea of subtracting a simple auxiliary function from the singular integrand 
and adding this contribution again. This auxiliary function has to be chosen in such 
a way that it cancels all singularities of the original integrand so that the phase-space 
integration of the difference can be performed numerically, even over the singular regions 
of the original integrand. In this difference the original matrix element can be evaluated 
without regulators for soft or collinear singularities. The auxiliary function has to be 
simple enough so that it can be integrated over the singular regions analytically with 
the help of regulators, when the subtracted contribution is added again. This singular 
analytical integration can be done once and for all in a process-independent way because 
of the general factorization properties of squared amplitudes in the singular regions. At 
NLO several subtraction variants have been proposed in the literature [4-8], some of which 
are quite general; at next-to-next-to-leading order subtraction formalisms are still under 
construction [9]. 

The dipole subtraction formalism certainly represents the most frequently used sub- 
traction technique in NLO calculations. It was first proposed within massless QCD by 
Catani and Seymour [5] and subsequently generalized to photon emission off massive 
fermions [6] 1 and to QCD with massive quarks [7,8]. Among the numerous applications 
of dipole subtraction, we merely mention the treatment of the electroweak corrections to 
e + e~ — > 4 fermions [11], which was the first complete treatment of a 2 — > 4 particle process 
at NLO. The formulation [5,7,8] of dipole subtraction for NLO QCD corrections assumes 
so-called infrared safety of observables, i.e. that all soft or collinear singularities cancel 
against their counterparts from the virtual corrections, either after parton-density redefi- 

1 The case of light fermions, where masses appear as regulators, has also been worked out in Ref. [10]. 
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nitions for initial-state singularities or due to the inclusiveness of event selection criteria 
in soft or collinear configurations for final-state singularities. In Ref. [6] the collinear sin- 
gularities from photon radiation are regularized by physical fermion masses, and only for 
final-state radiation these were assumed to cancel due to inclusiveness of the observable. 

In the following we generalize the method of Ref. [6] by dropping the latter assumption 
and by considering also other collinear-singular configurations involving photons, which 
are regularized by physical fermion masses: 

1. In Section 2 we deal with non-collinear-safe final-state radiation off light 
(anti-)fermions f, where collinear singularities arise from the splitting /* — > f'y. 
Here and in the following, asterisks indicate off-shell particles. By non-collinear- 
safe radiation we mean that a collinear fermion-photon system is not necessarily 
treated as one quasi-particle, which by contrast is the case in any collinear-safe ob- 
servable. In collinear-safe situations, which are usually enforced by photon recombi- 
nation or a jet algorithm, singularities from final-state radiation cancel according to 
the well-known Kinoshita-Lee-Nauenberg (KLN) theorem [12]. Non-collinear-safe 
final-state radiation off a fermion /, in general, leads to corrections oc alnm/ that 
are enhanced by a logarithm of a small fermion mass nif. 

2. Section 3 is devoted to processes with incoming photons and outgoing light fermions. 
Here the collinear-singular splitting is 7 — > //*, i.e. if an outgoing ( ant i-) fermion 
/ is allowed to be scattered into the direction of the incoming photon, the cross 
section receives an enhancement oc lnmj from this phase-space region. 

3. In Section 4 we treat processes with light fermion-antifermion pairs in the final state, 
i.e. when an outgoing photon with low virtuality splits into an // pair, 7* — > //. If 
the collinearly produced // pair can be distinguished from a plainly emitted photon 
(that has not split), the considered cross section again receives an enhancement 
oc mm/. 

4. Finally, in Section 5 we concentrate on processes with forward-scattered light 
(anti-)fermions, where the splitting / — > fj* leads to a collinear singularity if 
the emitted photon is almost real. Again this phase-space region enhances the cross 
section by a factor oc mm/. 

While Section 2 builds on the conventions and results of Ref. [6], Sections 3, 4, and 5 are 
self-contained and can be read independently. 

Of course, the considered situations could all be treated by fully including a non-zero 
fermion mass m/ in the calculation. However, if m/ is small compared to typical scales in 
the process, which is the case for electrons or muons in almost all present and future high- 
energy collider experiments, such a procedure is very inconvenient. The presence of very 
small or large scale ratios jeopardizes the numerical stability of phase-space integrations, 
and mass terms significantly slow down the evaluation of matrix elements. The subtraction 
technique described in the following avoids these problems by completely isolating all 
mass singularities from squared matrix elements, so that finally only amplitudes for a 
massless fermion / are needed. We support particle polarization whenever relevant, in 
particular for all incoming particles. In order to facilitate cross-checks in applications, the 
corresponding formulas for the phase-space slicing approach are also provided. 
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Figure 1: Generic diagrams for photonic final-state radiation off an emitter i with a 
spectator j or a in the final or initial state, respectively. 



In Section 6 we demonstrate the use and the performance of the methods presented in 
Sections 3, 4, and 5 in the example e~7 — > e~/i~/i + . A summary is given in Section 7, and 
the appendices provide more details on and generalizations of the formulas presented in 
the main text. In particular, the derivation of the factorization formulas for processes with 
incoming polarized photons splitting into light fermions and for the forward scattering of 
incoming polarized light fermions is described there. 



2 Non-collinear-safe photon radiation off final-state fermions 

2.1 Dipole subtraction and non-collinear-safe observables 

For any subtraction formalism the schematic form of the subtraction procedure to 
integrate the squared matrix element J2x \-Mi\ 2 (summed over photon polarizations A 7 ) 
for real photon radiation over the (N + l)-particle phase space d$i reads 

J d$x ]T \M!\ 2 = J d$i (j2 l-^il 2 - l-^subl 2 ) + /d$ ® (| [dA;] |.M sub | 2 ) , (2.1) 

where d$o is a phase-space element of the corresponding non-radiative process and [dk] 
includes the photonic phase space that leads to the soft and collinear singularities. The 
two contributions involving the subtraction function |.M su b| 2 have to cancel each other, 
however, they will be evaluated separately. The subtraction function is constructed in 
such a way that the difference ^a 7 |-Mi| 2 — |-M S ub| 2 can be safely integrated over d$i 
numerically and that the singular integration of |-M S ub| 2 over [dA;] can be carried out 
analytically, followed by a safe numerical integration over d$o- 

In the dipole subtraction formalism for photon radiation, the subtraction function is 
given by [6] 

|AW$i; K /)! 2 = - E Qf°fQf>°f'e 2 9ffi!(Pf,Pf', k) \M ($ ,//' 5 ™ f ) f , (2.2) 
/¥/' 

where the sum runs over all emitter-spectator pairs //', which are called dipoles. For a 
final-state emitter (final-state radiation), the two possible dipoles are illustrated in Fig. 1. 
The relative charges are denoted Qf, Qp, and the sign factors &f,crf> = ±1 correspond 
to the charge flow (<jf = +1 for incoming fermions and outgoing antifermions, aj = —1 
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for outgoing fermions and incoming antifermions). The implicitly assumed summation 

over r = ± accounts for a possible flip in the helicity of the emitter /, where Kf — ± is 

the sign of the helicity of / both in |.Mi| 2 and |.M su b| 2 . The singular behaviour of the 

subtraction function is contained in the radiator functions gf S J^l(pf,py, k), which depend 

on the emitter, spectator, and photon momenta pf, py, and k, respectively The squared 

lowest-order matrix element |.Mo| 2 of the corresponding non-radiative process enters the 

~(ff f ) ~(ff f ) 

subtraction function with modified emitter and spectator momenta pj and pf, . For 

(ff f ) (ff f ) 

a final-state emitter /, the momenta are related by pj + k ± py — pf ± py , where 
± refers to a spectator /' in the final or initial state, and the same set {k n } of remaining 
particle momenta enters |.Mi| 2 and |.Mo| 2 - The modified momenta are constructed in 
such a way that p K j — > pf + k in the collinear limit (pfk — > 0). 

Note that no collinear singularity exists for truly massive radiating particles /, because 
the invariant pyk does not tend to zero if the photon emission angle becomes small (for 
fixed photon energy k°). In such cases the corresponding masses are kept non-zero in 
all amplitudes, in the subtraction functions, and in the kinematics, and the subtraction 
procedure works without problems. Collinear (or mass) singularities result if the mass 
to/ of a radiating particle is much smaller than the typical scale in the process under 
consideration. In such cases it is desirable to set to/ to zero whenever possible. In a 
subtraction technique this means that to/ = can be consistently used in the integral 
/ d$i (j2x 1 \Mi\ 2 — |.M su b| 2 ), but that the readded contribution J[dk] |.M su b| 2 contains 
mass-singular terms of the form ahiTO/. If such mass singularities from collinear photon 
radiation do not completely cancel against their counterparts in the virtual corrections, 
the corresponding observable is not collinear safe. The dipole subtraction formalism as 
described in Ref. [6] is formulated to cover possible mass singularities from initial-state 
radiation, but assumes collinear safety w.r.t. final-state radiation. 

In collinear- safe observables (w.r.t. final-state radiation), and only those are considered 
for light fermions in Ref. [6], a collinear fermion-photon system is treated as one quasi- 
particle, i.e., in the limit where / and 7 become collinear only the sum pf + k enters 
the procedures of implementing phase-space selection cuts or of sorting an event into 
a histogram bin of a differential distribution. Technically this level of inclusiveness is 
reached by photon recombination, a procedure that assigns the photon to the nearest 
charged particle if it is close enough to it. Of course, different variants for such an 
algorithm are possible, similar to jet algorithms in QCD. The recombination guarantees 
that for each photon radiation cone around a charged particle / the energy fraction 

« = (2 - 3) 

is fully integrated over. According to the KLN theorem, no mass singularity connected 
with final-state radiation remains. Collinear safety facilitates the actual application of 
the subtraction procedure as indicated in Eq. (2.1). In this case the events resulting from 
the contributions of |-M su b| 2 can be consistently regarded as iV-particle final states of the 

non-radiative process with particle momenta as going into Ai ($0,//') , i- e - the emitter 

iff) ~(ff) (ff) 
and spectator momenta are given by pf , py , respectively. Owing to pf — > pf + k in 

the collinear limits, the difference J2x |-Mi| 2 — |.M su b| 2 can be integrated over all collinear 
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regions, because all events that differ only in the value of Zf enter cuts or histograms in 
the same way. The implicit full integration over all Zf in the collinear cones, on the other 
hand, implies that in the analytical integration of |.M S ub| 2 over [dk] the Zf integrations 
can be carried out over the whole Zf range. 

In non- collinear- safe observables (w.r.t. final-state radiation), not all photons within 
arbitrarily narrow collinear cones around outgoing charged particles are treated inclu- 
sively. For a fixed cone axis the integration over the corresponding variable Zf is con- 
strained by a phase-space cut or by the boundary of a histogram bin. Consequently, 
mass-singular contributions of the form alnrrif remain in the integral. Technically this 
means that the information on the variables Zf has to be exploited in the subtraction 
procedure of Eq. (2.1). The variables that take over the role of Zf in the individual dipole 
contributions in |.M S ub| 2 are called Zij and Zi a in Ref. [6], where / = % is a final-state emit- 
ter and j/aa final-/initial-state spectator. In the collinear limit they behave as z^ — > zi 
and Zi a z i- Thus, the integral f d$i (j2x 1 |-Mi| 2 — |-M su b| 2 ) can be performed over the 
whole phase space if the events associated with |.M S ub| 2 are treated as (N + l)-particle 
event with momenta pf — > Zfypj , pf — > pj, , and k — > (1 — Zffjpf . This can 
be formalized by introducing a step function cu t(p/, k,pf, {k n }) on the (N + l)-particle 
phase space which is 1 if the event passes the cuts and otherwise. The set {k n } simply 
contains the momenta of the remaining particles in the process. Making the dependence 
on 6 cut explicit, the first term on the r.h.s. of Eq. (2.1) reads 



/d*. 



Y,\Mi\ 2 ® C ut(Pf,k,p f >,{k n }) 



- E \M^, r \ 2 e^(z ff ^ f \(i-z fr )^ f \^ f \{k n }) 



(2.4) 



where we have decomposed the subtraction function |-M S ub| 2 into its subcontributions 
|-Msub,//'| 2 of specific emitter-spectator pairs //'. Apart from this refinement of the cut 
prescription in the subtraction part for non-collinear-safe observables, no modification in 
|.M su b| 2 is needed. Since its construction exactly proceeds as described in Sections 3 and 
4 of Ref. [6] , we do not repeat the individual steps in this paper. 

However, the modification of the cut procedure requires a generalization of the evalu- 
ation of the second subtraction term on the r.h.s. of Eq. (2.1), because now the integral 
over Zff implicitly contained in [dk] depends on the cuts that define the observable. In 
the following two sections we work out the form of the necessary modifications, where we 
set up the formalism in such a way that it reduces to the procedure described in Ref. [6] 
for a collinear-safe situation, while the non-collinear-safe case is covered upon including 
extra contributions. 

2.2 Final-state emitter and final-state spectator 

For a final-state emitter % and a final-state spectator j with masses mj and rrij the 
integral of g^y \pi,Pj, k) over [dk] is proportional to 

Gtr\P? 3 ) = <!//,( I' /A,) n^dz^g^ip^k), (2.5) 
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where the definitions of Sections 3.1 and 4.1 of Ref. [6] are used. There the results for 
Gf^\Pfj) with generic or light masses are given in Eqs. (4.10) and (3.7), respectively. In 
order to leave the integration over Zij open, the order of the two integrations has to be 
interchanged, and the integral solely taken over t/ij is needed. Therefore, we define 



G {suh) (P 2 z -) 



P 4 



2/2 (Zij) 



dy 



l 3 



(2.6) 



Note that no finite photon mass m 7 is needed in the function Qf^\pfp z) in practice, 



because the soft singularity appearing at z 
prescription in the variable z, 



(sub) , 



1 can be split off by employing a 



z) + \g^\p*,z) 



(2.7) 



This procedure shifts the soft singularity into the quantity G\^\p^), which is already 
known from Ref. [6]. Moreover, the generalization to non-collinear-safe integrals simply 



reduces to the extra term .... 
where the full ^-integration is carried out 



Gij U r\Pij, z ) i which cancels out for collinear-safe integrals 



For arbitrary values of m { and rrij a, compact analytical result of Qf^\P 2 -, z) cannot 
be achieved because of the complicated structure of the integration boundary. Note, 
however, that only the limit vrti — > of a light emitter is relevant, since for truly massive 
emitters no mass singularity results. The case of a massive spectator j is presented in 
App. A; here we restrict ourselves to the simpler but important special case m, = 0. 



In the limit m,i 
totically given by 



and mi 



the boundary of the y^ integration is asymp- 



:i 



P 2 z 



V2{z) = 1, 



and the functions and quantities relevant in the integrand gfj^ behave as 



Pik 



if. 

-f-yij, Rij(y) 



l-y, r ij (y) = l. 



The evaluation of Eq. (2.6) becomes very simple and yields 
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<l sub) 



(Pii,z) 



Pff(z) 
1-z, 



In 



Pt 

mf 



+(l + z) ln(l-z), 



where Pff(z) is the splitting function, 



(2.8) 



(2.9) 



(2.10) 



PffM 



1 + z 2 



(2.11) 



Equation (2.10) is correct up to terms suppressed by factors of m,. For completeness, we 
repeat the form of the full integral G < f^\P 2 j ) in the case of light masses, 



(sub) . 



7T 



p^) = c(p 2 , m 2 )-- + i, (r-V'r, 



i 

2' 



(2.12) 



6 



with the auxiliary function 

£( PW) = ln (|i) + - Iw^) + Img!), (2 ,3) 

which are taken from Eqs. (3.7) and (3.8) of Ref. [6]. 2 

Finally, we give the explicit form of the ij contribution \M. S ub,ij(^i)\ 2 to the phase- 
space integral of the subtraction function, 

J d$i \M suh ,ij(^i; Ki)\ 2 = -^QiOiQjOj J d$ ,ii J dz 

} (2-14) 

X \Mo(pi,pj;TKi)\ 2 @ C ut(pi = zpi,k = (1 - Z)pi,pj, {kn}), 

generalizing Eq. (3.6) of Ref. [6]. While pi,pj, {k n } are the momenta corresponding to the 
generated phase-space point in $o,ij, the momenta pi and k result from pi via a simple 
rescaling with the independently generated variable z. The invariant is calculated via 
Pfj = (j>i+pj) 2 independently of z. The arguments of the step function B C ut(Pi, k, pj, {k n }) 
indicate on which momenta phase-space cuts are imposed. 

For unpolarized fermions the results of this section have already been described 
in Ref. [13], where electroweak radiative corrections to the processes 77 — > WW — > 
4 fermions were calculated. In this calculation the results for non-collinear-safe differential 
cross sections were also cross-checked against results obtained with phase-space slicing. 
Another comparison between the described subtraction procedure and phase-space slic- 
ing has been performed in the calculation of electroweak corrections to the Higgs decay 
processes H — > WW/ZZ — > 4 fermions [14]. 

2.3 Final-state emitter and initial-state spectator 

For the treatment of a final-state emitter % and an initial-state spectator a, we con- 
sistently make use of the definitions of Sections 3.2 and 4.2 of Ref. [6]. In this paper we 
only consider light particles in the initial state, because the masses of incoming particles 
are much smaller than the scattering energies at almost all present and future colliders. 
Therefore, the spectator mass m a can be set to zero from the beginning, which simplifies 
the formulas considerably. 

Before we consider the non-collinear-safe situation, we briefly repeat the concept of 
the collinear-safe case described in Ref. [6]. Following Eqs. (4.24) and (4.27) from there, 
the inclusive integral of gfa^\pi-,Pa-, k) over [dk] is proportional to 

Gt?{Pl)= j^dxGf^iPlx) (2.15) 

with 

Gt?(Pl Xia) = /"^ dz m gtf(^Pa, k), (2.16) 

Z Jzi(x ia ) 

2 If dimensional regularization is used to regularize the soft singularity instead of a finite photon mass, 
the photon-mass logarithm in C has to be replaced according to ln(m^) — > (47r/i 2 ) c r(l + e)/e + 0(e), 
where D = 4 — 2e is the dimension and \i the reference mass of dimensional regularization. 
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x\G^(P*)5(l-z)+\g^>(P*,z) 



(sub) , 



where we could set the lower limit x of the ^-integration to zero because of m a = 
0. Since, however, the squared lowest-order matrix element |.Mo| 2 multiplying gf^ m 



Eq. (2.2) depends on the variable Xi a , the integration of |.M su b| 2 over x = 
employing a [...]+ prescription, 



is performed 



X\ 



zi{x ia ) 
zi(x ia ) 



52 ) 

ia/ 



X) + 



? (sub) 
7 ia,r 



+ 



(2.17) 



This integration, where the ellipses stand for ^-dependent functions such as the squared 
lowest-order matrix elements and flux factors, is usually done numerically. Since the soft 
and collinear singularities occur at x — > x\ = 1 — (9(m 7 ), the singular parts are entirely 
contained in Gf^f > {Pf a ) in Eq. (2.17), and the upper limit X\ could be replaced by 1 in the 
actual x-integration. For completeness we give the explicit form of the functions Gf^ 
and Gia^r^ m the limit rrii — > 0, 

J2 



G 



(sub) , 
io,+ 1 



a\PL\M 



71 

y + 1 - 



x 



2 In 



x 



X 



G 



(sub) 



Pi 



1 

2' 



(2.18) 



which are taken from Eqs. (3.19) and (3.20) of Ref. [6]. 

In a non-collinear-safe situation, the ellipses on the l.h.s. of Eq. (2.17) also involve 
^-dependent functions, as e.g. ^-functions for cuts or event selection. Thus, also the 
integration over z ia has to be performed numerically in this case, and we have to generalize 
Eq. (2.17) in an appropriate way. To this end, we generalize the usual [...]+ prescription 
in the following way. Writing 

/ d"r [g(r)]™ a) f(r) ee / d"rs(r) (/(r) - /(r)| ri=a ) (2.19) 

for the [...] + prescription in the reintegration in a multiple integral over n variables r k 
(k — 1, . . . , n), we can iterate this definition to two-dimensional integrals according to 



/ 



d"r 



(n,rj) 
+,(a,b) 



f(r) = / d"r 



in) 

+,(o) 



fa) 



,(&) 



= fd n rg(r)(f(r)-f(r) _ - f(r) + f(r 



(2.20) 



In the notation [<?(r)]+'( a ) we omit the superscript (r,) if g(r) depends only on the inte- 
gration variable r iy and we omit the subscripts (a) or (a, b) if a = 1 or a = b = 1. This 
obviously recovers the usual notation for the one-dimensional prescription used above. 
Introducing a double [. . .]+ prescription in x = x ia and z = z ia , we generalize Eq. (2.17) 
to 



XI 



zi{x) 



I JO Jz 1 (x) 



= £dx J o 1 dz{G^\Pl)5(l-x)5(l-z) + [g^(Plx)] + S(l-z 
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+ 



>(sub) 



7 ta,r V ta 



Pl,z) 5(1 -x) + 



,( sub ), 



x,z) 



(x,z) 



(2.21) 



If the functions hidden in the ellipses do not depend on z, the last two terms within the 
curly brackets do not contribute and the formula reduces to Eq. (2.17). 

We derive Eq. (2.21) and the explicit form of the two extra terms in two steps. In the 
derivation we quantify the previous ellipses by the regular test function /(x, z). The first 
step introduces a [. . .]+ prescription in the x-integration of the l.h.s. of Eq. (2.21) after 
interchanging the order of the integrations, 



I[f] = - 



P 2 

it 

2 Jo 



XI 



dx 



22 (z) 



pi i 

= — ^ / dz 

2 Jo 







la 

Y Jo 



dz 



zi(x) 
xi(z) 

xi (z) 



da^^a? f(x,z) 



dx 



(sub) 
9ia,T 



(x) 
+,(x 



i{z)) f(x,z)+g^ fMzlz)}. (2.22) 



The upper limit x\(z) of the x-integration follows upon solving the explicit form of the 
limits zi^(x) (given in Eq. (4.22) of Ref. [6]) for x. The full form of x\(z) is rather 
complicated for finite m 7 , but in the following it is only needed for m 7 = 0, where it 
simplifies to 



xi{z) 



(2.23) 



™ 7 =o Pf a z-m\{\-z)' 

Note that soft or collinear singularities result from the region of highest x values, x — > 
X\ = max{xi(z)}, so that the first term in curly brackets in Eq. (2.22) is free of such 
singularities owing to the [. . .] + regularization. Thus, we can set m 8 — > in this part, i.e. 
in particular x±(z) — > 1, yielding 



I[f\ = — f 



pi r i 



dz 



dx 



(sub) 
9ia,r 



(*) 



f(x, z) + f( Xl (z), z) j" (Z) dxgtf } • (2.24) 



In the second step we introduce a [...]+ prescription for the ^-integration in both terms, 



2 Jo 



dz 



dx 



(sub) 
9ia,T 



(sub) 
9ia,r 



(x) 
+ 



f(x, 1) 



pi -i 
ta / 

' 2 Jo 
pi -i 

ta / 

2 Jo 



+ f(xi(z),z) 

(x,z) 



rl 

f(x,z)+ / ( 
Jo 



dx / d^y^? 

JO L 

dz/(a;i(z),2) 



/ 

Jo 



ail(z) 



dx#) 



(sub) 



p? 

la 

~2~ Jo 

M 



(x) 



+ 



dx f(x, 1) 

P. 2 fl rxi(z) , 

/ dz / dxs£ 

2 Jo Jo 



sub) 



(2.25) 



In the second equality we just reordered some factors and integrations. Since all integrals 
over the test function / are now free of singularities, i.e. the singularities are contained 
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in the integrals multiplying /, we can set the regulator masses m 7 and m f to zero in the 
arguments of /. Thus, we can write 



9ia,r ( X , Z )\ + f{X,Z) + 



[ dxf(x,l) 
Jo 



+ 



(2.26) 



with the abbreviations 



-(sub)/ \ 
9ia,r ( X , Z ) 



Pi 



(sub) 



c (sub)/p2 
Y (sub) / p,2 



o 9ia,r 



2 jo 



P. 2 rxi(z) , . 

2 Jo 



GSS? ) (^) = -^/ o 1 dz / q 



xi(z) 



m^=0 



da; o^ sub) 



(2.27) 



Equation (2.26) is equivalent to the anticipated result (2.21), which was to be shown. 
The explicit results for Qf^\pf a) x) and Gf^(Pf a ) have already been given above in 
Eq. (2.18), the two remaining functions are easily evaluated to 



-(sub)/ \ 
9ia,+ ( X , Z ) = J 

Ot?(PLz) = i-z. 



x \2 — x — z 
Plz 



In 



- 1 - - ) 9i™-\ x , z ) = 0, 
2m(2-z) 



l-z 



+ (l + z) \n(l-z), 



(2.28) 



The collinear singularity oc lnm^ that appears in non-collinear-safe observables is con- 
tained in the function {Pia-> z ) ■ 

The resulting ia contribution |A^ S ub,ia( ( ^ ) i)| 2 to the phase-space integral of the subtrac- 
tion function reads 

Jd®! |-Msub,ia($i;^)| 2 = -^QaVaQiVi dx J d$ ,ia(P£, x) f dz 

x C ut(pi = zpi(x), k = (1 - z)pi(x), {k n (x)}) 



x - Gt™\pi) 5(1 -x) 5(1 -z)+ [Gt::\pL *)\ ai- Z ) 

X 1 J + 



+ 



QiZr\PLz)\5(l-x) + 



+ 



-(sub). 



(x,z) 



Mo(pi(x),p a (x);TK i )\ 2 , (2.29) 



which generalizes Eq. (3.18) of Ref. [6]. Again, the arguments of the step function 
@cut(Pij k, {k n }) indicate on which momenta phase-space cuts are imposed. We recall 
that <3>o,ia is the phase space of momenta pi(x) and {k n (x)} (without final-state radia- 
tion) with rescaled incoming momentum p a (x) = xp a instead of the original incoming 
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momentum p a . In the actual evaluation of Eq. (2.29), thus, the two phase-space points 
&o,ia(Pi a i x ) an d &o,ia(Pi a , x — 1) have to be generated for each value of x owing to the 
plus prescription in x. The relevant value of the invariant Pf a is then calculated sepa- 
rately via Pf a = (pi— pa) 2 for each of the two points, so that Pf a results from the momenta 
entering the matrix element A4q in both cases. 3 The variable z, however, is generated 
independently of the phase-space points and does not influence the kinematics in the 
matrix element. 

The combination of the subtraction procedures described in this and the previous 
section has been successfully applied and compared to results obtained with phase-space 
slicing in the calculations of electroweak corrections to Drell-Yan-like W-boson produc- 
tion, pp — > W — > v\l + X, and to deep-inelastic neutrino scattering, u^N — > u^/ ji + X, 
building on the calculations discussed in Refs. [15,16] and [17], respectively. 



2.4 Phase-space slicing 

In the phase-space slicing approach the soft and collinear phase-space regions are 
excluded in the (numerical) integration of the squared amplitude of the real-emission 
process. In the so-called two-cutoff slicing method the soft region is cut off by demanding 
that the photon energy k° should be larger than a lower cut AE which is much smaller than 
any relevant energy scale of the process. The collinear regions are excluded by demanding 
that each angle of the photon with any other direction of a light charged particle should 
be larger than the cut value AO <C 1. Note that this phase-space splitting is not Lorentz 
invariant. In the soft and collinear regions the photon phase space can be integrated out 
analytically by employing the general factorization properties of the squared amplitudes, 
which are, e.g., discussed in Section 2.2 of Ref. [6] (including polarization effects). General 
results for the integral over the soft region can, e.g., be found in Refs. [18,19]. The integrals 
over the collinear regions for final-state radiation can be easily obtained from intermediate 
results of the two previous sections as follows. 

The cuts defining the collinear region for the photon-emitter system of Section 2.2 
translate into new limits for the integration variables and Zij, 



1 - Zi 



P 2 z- 

1 13^13 



< Va < 



(ft ) 2 1 - z 



p?. 



A0 



< z^ < 1 - 



AE 



(2.30) 



which are asymptotically valid up to the relevant order in Wj — >• 0. With these new limits 
on i/ij we evaluate the integral defined in Eq. (2.6) and obtain 



S? t >{p°,z) = P„(z) 



V mf 



&f>(plz) = l-z. 



(2.31) 



The integrals of these functions over z = z i3 - are given by 

(AE 2 \ 3 



Gf>(p?) = - 



In 



2 

mf 



(2.32) 



3 For a more formal explanation of this subtle but important point we refer to the discussion at the 
end of Section 6.3 of Ref. [8]. 



11 



As it should be, in these results the dependence on the spectator particle j completely 
disappears, because it was only needed in the phase-space parametrization. We also note 
that the same results can be obtained from Section 2.3, where the limits on x ia and z ia 
are changed to 



, <l-x ia < ^-j ^ A0 2 , < ^ < 1 - —. (2.33) 

*ia Z 'ia~l~ m i{l- Z ia) z ia Pi 

Using the functions and G^, the integral over the collinear photon emission 
cone around particle % reads 

jf^.d^lMx^!;^)! 2 = |^Q 2 / d$ jTdz {<7W(p°)*(l-*) + [<^ sIi) (P^)] + } 

x |Alofe; ™;)| 2 ©cut(ft = zpi, k = (1 - z)pi, {&„}), (2.34) 

where the momenta and {/c n } belong to the phase-space point $ - Of course, apart from 
the polarization issue this is a well-known result which can be found in various papers [3]. 4 

3 Collinear singularities from 7 — > ff* splittings 

3.1 Asymptotics in the collinear limit 

We consider a generic scattering process 

7 (M 7 ) + a( Pa )^/( P/ )+X, (3.1) 

where the momenta of the particles are indicated in parentheses and A 7 = ± is the photon 
helicity. Here a is any massless incoming particle and / is an outgoing light fermion or 
antifermion. The remainder X may contain additional light fermions which can be treated 
in the same way as /. For later use, we define the squared centre-of-mass energy s, 

s = {p a + kf = 2 Pa k. (3.2) 

The collinear singularity in the squared matrix element \A4 ia ^fx\ 2 occurs if the angle Of 
between / and the incoming 7 becomes small; in this limit the scalar product (kpf) is 
of (9(m 2 ), where to/ is the small mass of /. Neglecting terms that are irrelevant in the 
limit to / — > the squared matrix element \Ai-y a ^fx{k,p a ,Pf] A 7 )| 2 for a definite photon 
helicity A 7 = ± (but summed over the polarizations of /) asymptotically behaves like 

\M ja ^fx(k,p a ,pf, A 7 )| 2 fe ~ Qje 2 h1 f (k,p f ) \Mf a ^ x (Pf = xk,p a ; kj = rA 7 )| 2 , (3.3) 

where x = 1 — p°j/k° and Qje is the electric charge of /. The matrix element M.j a ^ x 
corresponds to the related process fa^X that results from 7a(— * ff*a) — > fX upon 
cutting the /* line in all diagrams involving the splitting 7 — > //* (see also Fig. 2). The 
incoming momenta relevant in the different matrix elements are given in parentheses. 



descriptions of phase-space slicing for initial-state radiation off unpolarized particles can also be 
found in Ref. [3]; the case of polarized incoming particles is, e.g., treated in Ref. [20]. 
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Pj 



Figure 2: Generic diagrams for the splittings 7 — > //* with an initial-state spectator a or 
a final-state spectator j. 

Moreover, in Eq. (3.3) we assume a summation over r = ±, where r = ± refers to the two 
cases where the sign nj of the / helicity is equal or opposite to the photon helicity A 7 . 
The functions h^(k,pf), which rule the structure of the collinear singularity, are given 
by 

with the splitting function 

p h (x) = (1 - xf + x 2 . (3.5) 

The derivation of this result is given in App. B.l. 

Note that the collinear singularity for kpf — > can be attributed to a single external 
leg (namely /) of the related hard process fa—>X. Thus, there is no need to construct 
the subtraction function |-M S ub| 2 from several dipole contributions oc QfQf. Instead we 
can construct |.M su b| 2 as a single term oc Q 2 . Nevertheless we select a spectator /' to the 
emitter / for the phase-space construction, which proceeds in complete analogy to the 
photon radiation case. We have the freedom to choose any particle in the initial or final 
state as spectator. In the following we describe the "dipole" formalism in two variants: 
one with a spectator from the initial state, another with a spectator from the final state. 
The two situations are illustrated in Fig. 2. 

3.2 Initial-state spectator 

The function that is subtracted from the integrand \M.~, a -^fx{k,Pa,Pf] A T ) | 2 is defined 
as follows, 

|A^sub(A 7 )| 2 = Q 2 f e 2 h^ a (k, Pf , Pa ) \M f - a ^ x (p f -,Pa,{kn};Kf- = T\ 7 )\ 2 , (3.6) 
with the radiator functions 

hf a (k, Pf , Pa ) = — i_ (P fl (x f , ia ) + X -^l) - h^(k,p f , Pa ), 



Xf,ya(kp f ) V kp f 
hl> a (k,p hPa ) = -^^ {1 - Xffia ){l- Xftia - ), (3.7) 
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and the auxiliary quantity 

p a k - p f k - p aPf 
/l7a ~ • ^ ' 

Pa^ 

Here we kept the dependence on a finite rrif, because it is needed in the analytical treat- 
ment of the singular phase-space integration below. The modified momenta pj and {k n } 
entering the squared matrix element on the r.h.s. of Eq. (3.6) will only be needed for 
rrif = in applications with small values of rrif. In this limit they can be chosen as 

p f j(x)=xk^ py = pJ(x />7B ), K = WvK ( 3 - 9 ) 
with the Lorentz transformation matrix A ll v given by 

Aw,- (p+ ^ + ( y>» + ^. (»o) 

= p£ + F - pj, P»(x)=pZ+p^x), F»=rt+pr}. (3.11) 

It is straightforward to check that |.M su b| 2 possesses the same asymptotic behaviour 
as \Ai ia ->fx\ 2 in Eq. (3.3) in the collinear limit with rrif — * 0. Thus, the difference 
|.M 7 a->/x| 2 — |-M su b| 2 can be integrated numerically for rrif = 0. 

The correct dependence of |-M S ub| 2 (and the related kinematics) on a finite rrif is, 
however, needed when this function is integrated over Of leading to the collinear singularity 
for 9f — > 0. The actual analytical integration can be done as described in Ref. [6] (even 
for finite m a and rrif). Here we only sketch the individual steps and give the final result. 
The (N + l)-particle phase space is first split into the corresponding iV-particle phase 
space and the integral over the remaining degrees of freedom that contain the singularity, 

J d(f)(pf,P;k + p a ) = j dx J d(j>(P(x);pf(x)+Pa) J [dp f (s,x,y f ^ a )}, (3.12) 

with the explicit parametrization 

/s rw2( x ) r 
[dp f {s,x,y Lja )} = -^—^ dy fr/a J d<f> f . (3.13) 

The upper kinematical limit of the parameter x = #/ j7a is given by 

*i = l-^, (3.14) 

V s 

but in the limit rrif — > we can set x\ = 1. While the integration of the azimuthal angle 
(f)f of / simply yields a factor 2n, the integration over the auxiliary parameter 

„ _ k Pf _ 2k Pf fQ 1 

V ^-kp- a -— (3 - 15) 



with the boundary 

yiA x ) = \ 



. . 4m 2 f 
1 - x =F \ (1 - x) 2 - 



(3.16) 
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is less trivial. Defining 

ry 2 (x) 
lyi(x) 

the result of this straightforward integration (for mj — > 0) is 



W/> a (s,x) = ^ f MX) dy f , ia h^ a (k,p f ,p a ), (3.17) 

Z fa;) 



Hl f ' a (s,x) = P fy (x) ln( s(1 ^/ )2 ^) + 2x(l -x)- Hl f ' a (s, x) 



Hl f ' a (s,x) = (l-x) 2 ln^ S[1 m2 ^ -(l-x) 2 . (3.18) 

For clarity we finally give the contribution a s ^i^fx that has to be added to the result for 
the cross section obtained from the integral of the difference \M ia ^fx\ 2 — |-M su b| 2 , 



a s ^ fx (k,p a ; A 7 ) = N cJ jf dxH^ a (s,x) a fa ^ x (p f = xk,p a ; Kf = rA 7 ). (3.19) 

Although formulated for integrated cross sections, the previous formula can be used to 
calculate any differential cross section after obvious modifications. 

For the case of unpolarized photons this subtraction variant has already been briefly 
described in Ref. [17], where it was applied to the contributions to deep- inelastic neutrino 
scattering, u^N — > v^j ' ji + X, that are induced by a photon distribution function of the 
nucleon N. Moreover, the method presented here was successfully used in the calcula- 
tion of photon-induced real corrections to Drell-Yan-like W production (see Section 10 
of Ref. [1] and Ref. [16]) and of photon- and gluon-induced real corrections to Higgs pro- 
duction via vector-boson fusion at the LHC [21]. All these results were also cross-checked 
against phase-space slicing. 

3.3 Final-state spectator 

As an alternative to the case of an initial-state spectator described in the previous 
section, we here present the treatment with a possibly massive final-state spectator j 
with mass rrij, i.e. we consider the process 

7 (fc) + a( Pa ) - f(p f )+j( Pj )+X. (3.20) 

The initial-state particle a is assumed massless in the following, but all formulas can 
be generalized to m a ^ following closely the treatment of phase space described in 
Section 4.2 of Ref. [6]. The subtraction function now is constructed as follows, 

|M sub (A 7 )| 2 = Q)e 2 h] f T (k,p f ,p 3 ) \Mf^j X {pf,Pa,Pj\ Kf = r\j)\ 2 , (3.21) 
with the radiator functions 

h&(k, Pf , Pi ) = — -|—y (P»(x tj „) + X -^) - ht(k,p f , Pj ), 



xfiAkPf) V /7V Ihl> kp f 

hl(k, Pf ,p 3 ) = -* fl- Xfj „)(l- Xfj „-^ (3.22) 
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and the auxiliary parameter 



x fjn 



k Pj + kp f - PfPj 



(3.23) 



kpj + kpf 

The momenta pj and pj are given by 

f f {x)=xk», pj- = pj-(a; /jl7 ), g = P"+pJ-, P»=rf+p$-hf, (3.24) 

while the momenta of the other particles are unaffected. Note that this construction of 
momenta is based on the restriction rrif — 0, which is used in the integration of the 
difference |A4 7a ^/jx| 2 — |A^ su b| 2 - 

In the integration of |-A4 S ub| 2 over the collinear-singular phase space, of course, the 
correct dependence on a finite to/ is required. Owing to the finite spectator mass rrij, this 
procedure is quite involved; we sketch it in App. B.2. Here we only present the results 
needed in practice. The cross-section contribution a^^/jx that nas to be added to the 
integrated difference \M.-, a ->fx\ 2 — \-M S ub\ 2 is given by 

^ X tfjx(k, Pa ;\,)=N cJ ^ j\xH] f T {P\x)o Ia ^ x {p f = xk,p a ;K f - = rA 7 ), (3.25) 
where the collinear singularity is again contained in the kernels 

H] f + (P 2 ,x) = -P f ,(x)ln 



mjx 



m 2 x 



m 2 - P 2 )(l-x) V (to 2 -P 2 ){1 -x) 



n]Up\x), 



H]f_(P 2 ,x) = -(l-x) 2 \n 



mjx 



to: 



P 2 )(l-x) 



1 + 



m 2 x 



TO' 



P 2 )(l-x) 



+ 2x(l - x) 

-(1-x) 2 . 

(3.26) 



Of course, the singular contributions oc In toj have the same form as in the case of an 
initial-state spectator discussed in the previous section. 

3.4 Phase-space slicing 

From the results of the two previous sections, the corresponding formulas for the phase- 
space slicing approach can be easily obtained. The collinear region, which is omitted in 
the phase-space integration, is defined by the restriction Of < A6 on the fermion emission 
angle Of in some given reference frame. 

In Section 3.2 this constraint translates into new limits on the variable y/, 7a , 



m, 



(k ) (I Xjyra) Aa 2 

< Vf,ja < ; ^ , 



S(l - Xf„ a ) S 

which modifies the result of the integral analogously defined to Eq. (3.17) to 

•(A; ) 2 (l-:r) 2 A^ 



(3.27) 



Hl f (k°,x) = P fl (x)\n 



TO 



+ 2x(l - x) - Hl f (k°, x), 



Hl f {k°,x) = (f-a:) 2 ln 



/ 

-x) 2 A6' 



TO 



-(1-xf 



f 
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(3.28) 



The cross-section contribution of the collinear region of / then reads 

d$Hfx%Pa\ A 7 ) = N cJ J dxH1 f (k°, x) <Jj a ^ x {Pj = xk,p a ; kj = rA 7 ). (3.29) 

The same result is obtained from Section 3.3 with App. B.2, where the new limits on 
zjj^ read 



4 Collinear singularities from 7* — ► // splittings 

4.1 Asymptotics in the collinear limit 

We consider a generic scattering process 

a( Pa ) + b( Pb ) -> + /(p 7 ) + X, (4.1) 

where the momenta of the particles are indicated in parentheses. Depending on the par- 
ticle content of the remainder X, there may be additional, independent collinear-singular 
configurations, but we are interested in the region where the invariant mass (pf +Pf) 2 = 
2m 2 + 2pfPf of the produced fermion-antifermion pair // becomes of the order 0(m 2 ), 
where m/ is small compared to typical scales in the process. The singular behaviour 
of the full squared matrix element \-M a b-,ffx(.PfiPf)\ 2 entirely originates from diagrams 
containing a 7* — > // splitting, i.e. the singularity is related to the subprocess ab — > 7X. 
For the matrix element of this subprocess we write Ai a b^x = T^ b ^ /X (k)e\ Jifl (k)* , where 
Tj^ x {ti) is the amplitude without the photon polarization vector S\ itiJi {k)*. In the 
collinear limit PfPj — > the light-like momentum k is equal to k = p/ + pj up to mass- 
suppressed terms. Neglecting terms that are irrelevant in the limit m/ — > the squared 
matrix element \M. a b^ffx(PfiPf)\ 2 asymptotically behaves like 



\M ab ^ ffx (Pf,Pf)\ 2 P/ ~ N c jQje 2 h ff -^(p f ,p f )T^ jX (krT^ lX (k), (4.2) 



where 



h f f,^(p f ,Pf) 



(Pf+Pf) 2 [ ^ — m^J ' A; ' 



* = % (4-3) 



and N c j is the colour multiplicity of / (iV Cj i epton = 1, A c quar k = 3). The momentum k± 
is the component of P f that is orthogonal to the collinear axis defined by k, i.e. kk± = 0, 
and becomes of 0{mf) in the collinear limit. An explicit prescription for the construction 
of k± can, e.g., be found in Ref. [8], where the analogous case of the gluonic splitting 
into massive quarks Q, g* — > QQ, is worked out. It is important to realize that hfj^ 
in Eq. (4.2) is not proportional to the polarization sum E^ v = ^x^,fj,(k)*S\^,u(k) of 
the photon, so that the r.h.s. is not proportional to the polarization-summed squared 
amplitude \M. a b^- 1 x\ 2 of the subprocess. This spin correlation has to be taken care of 
in the construction of an appropriate subtraction function in order to guarantee a point- 
wise cancellation of the singular behaviour in the collinear phase-space region. The spin 
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Figure 3: Generic diagrams for the splittings 7* — > ff with an initial-state spectator a or 
a final-state spectator j, where / is a light fermion or antifermion. 

correlation encoded in htf „ v drops out if the average over the azimuthal angle <pf of the 
7* — > ff splitting plane around the collinear axis is taken. 5 Indicating this averaging by 
(. . .) ( f )f = f d0//(27r), we get {hfj^ u )^ f = E^ u hfj with (in four space-time dimensions) 



2 . . 2m 2 , 



(4.4) 



up to terms that are further suppressed by factors of ntf. The averaged squared matrix 
element behaves as 



{\M ab _ ff - x (pf,p f -)\ 2 }^ f 



Pfpp 



N cJ Qje 2 h fI (p f ,p f -) \M ab ^ x (k)\ 2 . (4.5) 



Since the collinear singularity for PfPj — > can be attributed to a single external leg 
(the photon) of the related hard process ab —>■ 7X, also in this case there is no need to 
construct the subtraction function | -M su b| 2 from several dipole contributions. The function 
|-M S ub| 2 can be chosen as a single term oc Q 2 . Nevertheless a spectator is selected for the 
phase-space construction, as in the previous section. In the following we describe the 
"dipole" construction in two variants: one with a spectator from the initial state, another 
with a spectator from the final state. The two situations are illustrated in Fig. 3. 

4.2 Initial-state spectator 

We define the subtraction function as 



l-Msubl 2 = N c jQ 2 f e 2 hy^(p f} p f -,p a )T^ x (p a ,kyT^^ 



(4.6) 



with 



tiff(Pf,Pf,Pa) 



(pf + PfY 



/J.V 



'Pf ~ Z ff,aPf 

(4.7) 



5 As described in Refs. [5,8], for unpolarized situations this average can be easily obtained upon con- 
traction with the projector ^d^k) = ^[— i g Al ^ + ( "gauge terms" involving fc M or fc I/ )]/(l — e), which fulfills 
—g^ v d^ v (k) = 2(1 — e) and fc M <i M „(fc) = in D = 4 — 2e space-time dimensions. 
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and the auxiliary parameters 



x 



ff,a 



PaPf + PaPf ~ PfPf ~ mj 
PaPf + PaPf 



C ff,a 



l-Z 



ff,a 



PaPf 



PaPf + PaPf 



(4.8) 



The auxiliary momenta entering the amplitude for the related process ab — > 7X are given 
by 

~k»(x) = P» + %(x), P = fc"(s //>0 ), P» = p» } + p) - pS, (4.9) 

while the momenta of the other particles remain unchanged. In these equations we kept 
the dependence on rrif, but of course in the numerical integration of \M. a b^fjx\ 2 ~ |-Msub| 2 
we can set to/ to zero, because we are only interested in the limit rrif — > 0. For the 
integration of |-M su b| 2 over the collinear-singular region, we need the TO/-dependence of 
the spin average of h a ^ u , 



2m) 



(4.10) 



h° f -( Pf , Pf -, Pa ) = \P h {z fLa ) + 

and an appropriate phase-space splitting, 

J d<f>(p f , Pr ,P + p a ) = J Q Xl dx J dc!>(k(x);P + p a (x)) J[dp f (P 2 ,x,z)], (4.11) 

where we have used the shorthands x — Xjj a and z = Zff a . The explicit form of J[df>/] 
reads 

r n —P 2 1 rMx) r 

with the integration limits for the variables x and z 

p2 



Xi 



p 2 -4 m y 



1 



1 ± 



/ Xi — x 
Xi(l — x) 



(4.13) 



Separating the singular contributions as described in Section 2.3, we rewrite the integral 
of h'j j for rrif — > as 



/ da; 

2 Jo 



Z2{x) 

zi(x) 

1 



dz h a f f-(pf,p f -,p a )--- 



dx I dz \ H a ff (P 2 ) 5(1 - x) 5(1 -z)+ H a ff (P\ x) : 5(1 - z) 



+ 



n a ff (P 2 ,z)\ 5(l-x) + [h a ff -(x,z) 



(x,z) 

+ 



(4.14) 



The new functions h a ^, etc., defined here are obtained from obvious substitutions and 
straightforward integrations, 



h} f (x,z) = - [ ^Pf 1 (z) 
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n W p2 i x ) - 3(1 _ x y 



n a fI (p\z) = p h (z) 



In(- P2 ^)-1 



m) 



+ 2z(l-z), 



-P 2 
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(4.15) 



Using these functions the phase-space integral of the subtraction function reads 



&® f! \M snh (® f j)\ 2 = N cJ ^£- £dx |d$ 
x 6 cu t(p/ = zk(x),pf = (1 - z)^(x), 



Jo 



x - \ H a ff -(P 2 ) 5(1 - x) 5(1 -z)+ [n a ff -(P\ x) 5(1 - z) 

X I 



+ 



H a s] (P\z)\ 5(1 -x) + [h a ff (x,z)] {X ' Z) \ \M ab ^x(pa(x),k(x))\ Z , (4.16) 



where we have made explicit which momenta enter the cut function Q C ut(Pf,Pf, {k n })- 
Concerning the phase-space integration over d<l ) 7 (P 2 , x) and its integration over the boost 
parameter x the same comments as made after Eq. (2.29) apply. There are actually two 
phase-space points for each x value to be generated (one for x < 1 and another for 
x = 1), each determining momenta p a (x), k(x), {k n (x)} for the evaluation of P 2 and the 
matrix elements. The generation of the parameter z proceeds independently, and the 
squared amplitude \M. a b~*-yx\ 2 in Eq- (4.16) does not depend on z. Thus, if the full range 
in z is integrated over, i.e. if the collinear // pair is treated as a single quasiparticle 
in the cut procedure, the last two terms in curly brackets do not contribute. In this 
case the fermion-mass logarithm is entirely contained in the H^j contribution. According 
to the KLN theorem this contribution will be completely compensated by virtual 0(a) 
corrections to the process ab — > 7X if collinear // pairs are not distinguished from emitted 
photons. 

4.3 Final-state spectator 

Since the case with a massive final-state spectator j is quite involved, we here present 
the formalism for rrij = and give the details for the massive case in App. C. 
For rrif = rrij = 0, the subtraction function can be defined as 



l-Msubl = N cJ Q f e h ff -^ flu (p f ,p f -,p j )T^ ljX (k,p j yT^ ljX (k,p j ) 



with 



and the auxiliary parameters 



Pfpf 



i'-f fjl''r - ZffjPf) {zffjP 1 } ~ 1 



*ffjPf 



z ffj ~ 1 ~ z ffj ~ 



PfPj 



PfPj + PfPj 



Vffj = 



PfPf 



PfP 3 + PjP 3 + PfPf 



(4.17) 



(4.18) 



(4.19) 
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The new momenta entering the amplitude for the related process ab — * 7j'X are given by 



Pi 



Pj/i 1 - Vff. 



3>i 



(4.20) 



whereas all remaining momenta k n of particles in X remain unchanged. Equation (4.17) 
can be used to integrate the difference \M. ab ^fjj X \ 2 — |-M su b| 2 fc> r massless fermions /. In 
order to integrate |-M su b| 2 over the collinear-singular region, the dependence on rrif has 
to be taken into account. Details of this procedure can be found in App. C. The result 
can be written in the form 



d$ // |A^ sub ($ // )| 2 = N cJ ^ J d$ 7 £dze cut ( Pf = zk,Pf 



1 - z)k,pj, {k n } 
\M ' ~- " l2 



(4.21) 



with 



n m {p\z) = p Sl {z) 



In 



P 2 z(l - z) 



m, 



- 1 



+ 2z(l-z), 



HffAP 2 ) 



Km) 
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9 



(4.22) 



The momenta k,pj,{k n } directly correspond to the generated phase-space point in $ 7 , 
while the parameter z is generated independently. The comments on the ^-integration 
made at the end of the previous subsection apply also here. The squared amplitude 
\M-ab^~jjx\ 2 in Eq. (4.21) does not depend on z, and thus, if the event selection for / and 
/ is inclusive in the collinear region of the 7* — > // splitting, the integral over z trivially 
reduces to the factor Hfjj(P 2 ). 

4.4 Phase-space slicing 

Here we again deduce the integral over the collinear phase-space region which is needed 
in the slicing approach. This region can, e.g., be defined by restricting the angle Qfj 
between the / and / directions to small values, 6fj < A0 <C 1. 

In Section 4.2 this restriction leads to new limits in Xfj a and Zff a , 



m. 



(k 



0\2 



-P 2z ff,a{ 1 ~ Z ff,a) 



< 1 ~~ X ff,a < ^ZpZ'ff*- 



ZffM - z ff - a )A0 2 , < z ff ~ < 1, (4.23) 



where k° = p° + p°j is the energy in the // system. This modifies the integrated results 
to 



n ff -{k°,z) = P f7 (z)\n 



(k°) 2 z 2 (l - z) 2 A6> 5 



ml 



+ 2z(l-z), 



2 ( (k°) 2 A9 2 \ 23 



(4.24) 



where TLfj and Hfj are defined analogously to Eq. (4.14). The integral of the squared 
matrix element over the collinear regions then reads 
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I 



x e cut (p/ = zk,p f = (1 - z)k, {k n }) (4.25) 
x {H ff -(k°)5(l-z)+ [n fI (k°,z)}^\M ab ^ x (k)\ 2 . 

The same results can be obtained from Section 4.3 with App. C, where the new limits 
on the integration variables are given by 



i' 2 W -*//;) ffJ P 2 



< Vffi < -nTZffli 1 ~ Wi)^, < *ffi < !■ ( 4 - 26 ) 



5 Collinear singularities from / — > fj* splittings 

5.1 Asymptotics in the collinear limit 

We consider a generic scattering process 

f(p f , Kf ) + a(p a )^f(p' f )+X, (5.1) 

with the momenta of the particles and the (sign of the) helicity Kf = ± of the in- 
coming fermion / indicated in parentheses. We are interested in the region where the 
squared momentum transfer {pf — p'^) 2 = 2m 2 — 2pfp'f of the scattered fermion / be- 
comes of the order 0(mj), where m/ is small compared to typical scales in the pro- 
cess. The singular behaviour of the full squared matrix element \M.f a -^fx(Pf,Pf', K /)| 2 
entirely originates from diagrams containing an / — > f'y* splitting, i.e. the singularity is 
related to the subprocess 7a — > X. For the matrix element of this subprocess we write 
M ya ^x{k,p a , A 7 ) = T^ x (k)e x ^^{k), where T^ x (k) is the amplitude without the pho- 
ton polarization vector e\^ tfM (k). In the collinear limit Pfp'f — > the momentum k is given 
by k = pf — p'f up to mass-suppressed terms. Neglecting terms that are irrelevant in the 
limit rrif — > the squared matrix element \M.f a ->fx(pf,Pf', K /)| 2 asymptotically behaves 
like 

\Mf a ^fx(Pf,Pa,p' f ; «/)| 2 pf ^ Nc >f Q)^ hf 4,Hu(Pf,P'f) T^ x (k,p a )* T^ x (k,p a ), 

(5.2) 

where 



4(1 - x) k ±jli k ±jV 
x 2 kj_ - x 2 mj 



x V (Pf-P'f) 

with 



+ V f 2 " X + 7^~%) {z+AkTe+Ak) ~ e-Ak)*e^(k 



k° , 
x = -0 . (5.4) 

Pf 
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Figure 4: Generic diagrams for the splittings / — > fj* with an initial-state spectator a or 
a final-state spectator j, where / is a light fermion or antifermion. 



The momentum k± is the component of k that is orthogonal to the collinear axis defined 
by pf, i.e. kxPf = 0, and becomes of 0(rrif) in the collinear limit. A derivation of this 
factorization is described in App. D.l. Note that h*J ^ in Eq. (5.2) is not proportional to 

the polarization sum E^ v = J2x £ x 1 ,fi(k)*ex ltU (k) of the photon, so that the r.h.s. is not 
proportional to the polarization-summed squared amplitude \M. ia _>x\ 2 of the subprocess. 
This spin correlation has to be taken into account in the construction of an appropriate 
subtraction function in order to guarantee a point-wise cancellation of the singular be- 
haviour in the collinear phase-space region. The spin correlation encoded in h{/ ^ v drops 
out if the average over the azimuthal angle (p'j- of the / — ► f'y splitting plane around the 
collinear axis is taken. Details of this averaging process, which is indicated by (. . .)^> , are 
given in App. D.l. The result is 

(l^/o-/x(p/,Pa,P/;«/)| 2 )^, p ^r_ tQ N cJ Q)e 2 h f T f (p f ,p' f ) \M ia ^x(k,p a ; A 7 = TK f )\ 2 

(5.5) 

with summation over r = ± and 



H f ip f ,p' f ) 1 



, x 2xm 2 f ( 2x 2 m 2 f 

P lf {x) + - f — + r[2-x+ f 



x(Pf-p' f ) 2 f 7JV ~' ' (pf-p' f ) 2 V" ~ (Pf-P'f) 2 ' 

which is valid in four space-time dimensions up to terms that are further suppressed by 
factors of m/. Here Pjf(x) is the splitting function 

P 7/ (x) = l + {l ~ X? - (5-7) 

Since the collinear singularity for p/p'f — * can be attributed to a single leg (the 
photon) of the related hard process 7a — > X, also in this case there is no need to construct 
the subtraction function |-M S ub| 2 from several dipole contributions. The function |A4 su b| 2 
can be chosen as a single term oc Q 2 , and a spectator is only used in the phase-space 
construction as previously. In the following we again describe the "dipole" construction 
in two variants: one with a spectator from the initial state, another with a spectator from 
the final state. The two situations are illustrated in Fig. 4. 

5.2 Initial-state spectator 

We define the subtraction function as 

\M suh (n f )\ 2 = N c jQ)e 2 h^ u {p }) p^p a )T^ x {k,p a y (5.8) 
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with 

h f K f ; a ^(Pf,p' f ,Pa) = 



-1 



4(l-x />/a ) l\k\ 



X fJa 



K _L X fja m f 



(P.f-P'f) 2 [ * 

+ K i{ 2 ~ x+ j^^) - 

and the auxiliary parameters 



(5.9) 



PaPf - PfP'f - PaP'f + m) 



Vfja = 



PfPf 



m 



PaPf PaPf 

Assuming again the incoming particle a to be massless and defining 

s = (p f + Pa) 2 =m} + 2 PaPf = s + m 2 , P» = jf f + p£ - pf, 
the needed auxiliary momenta for the related process 7a — > X are given by 



(5.10) 



(5.11) 



= a; (pi} 



in 



Pah 



fr = kr(x ftfa ), 



K _L — Pf 7 Pai 



(5.12) 



s --/■ - ,,„k 

P*(x) = P* = P*(x ftfa ), 

Z.H _ AM 

where the Lorentz transformation matrix A.^ v is constructed from the momenta P M and 
P M as in Eq. (3.10). In these equations we kept the dependence onmj, but of course in 
the numerical integration of \A4 f a ^fx\ 2 — |-M su b| 2 we can set rrif to zero if we are only 
interested in the limit rrif — > 0. For the integration of |-M S ub| 2 over the collinear-singular 
region, we need the m/-dependence of its azimuthal average, 

(|AW«/)| 2 )^ = N cJ Q 2 f e 2 h f T f ' a (p f: p' f:Pa ) \M ia ^x{k,Pa]\ = ™/)| 2 (5.13) 
with summation over r = ± and 

P lf {x) 



h f T f ' a (p f ,p' f ,p a ) 



sxy 



2x(l — x)m 2 



+ T12-X 



2x 2 m 2 f 



y[s(l - x - y) - m){2x + y)\ \ sy 

(5.14) 

where we have used the shorthands x = Xfj a and y = yfj a - An appropriate phase-space 
splitting is given by 



J d<P{Pf,P;Pf +Pa) = J dx J d<j)(P(x); k(x) +p a ) J [dp' f (s,x,y)} 
with the explicit form of / [dp^] 

/s ry?.( x ) r 
W 8 > x >vK = M2^l (x) dy J d(f) 'f 



and the integration limits for the variables x and y 



Xi 



y/s + mf' 



s l 2m 2 f 4m 2 f 
yi,2{x) = — 1 - X j^X =F \/ (1 - x) 2 - 



X 



(5.15) 



(5.16) 



(5.17) 
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In the limit mj — > the integral 

rr XS ryi(. X ) rr 

U^\s,x) = - dyh^ a (p f ,p' f , Pa ) 

Z Jyi(x) 

can be easily evaluated to 

■Vs(l-x) 



H f T f ' a (s,x) = In 



xrrif 



P, f (x) 



~(2-x) 



1-x 



X 



-r(l-x), 



and the part to be added to the cross section reads 

Q 2 f a 



v S ft^fx{PfiPa\K f ) = J dxH f T f ' a (s,x)a ia ^ x (k = xp f ,p a -\ = tkj). 



(5.18) 



(5.19) 



(5.20) 



5.3 Final-state spectator 

As an alternative to the case of an initial-state spectator, we now present the treatment 
with a final-state spectator j, i.e. we consider the process 

f(p ft K f ) + a( Pa ) -> /(//,) • j( Pj ) • .V. (5.21) 

The particles a and j are assumed massless in the following; the case of a massive spectator 
j is described in App. D.2. The subtraction function is constructed as follows, 

\M suh (K f )\ 2 = ^/^e 2 /i^(p / y / ,p J 0^^(*,Pa,Pi)*^ i j C (*,Pa,Pi) (5.22) 



with 



h 



-1 



2x 2 m 2 



( Pf - p' f yx 2 fj jz fj j 



(p/ -p'/) 2 



P/P/ 



P/P/ + P/Pj 



+ ^(2-a; + 

and the auxiliary parameters 

x = PfP'f + PfPj - P'fPj 
fKf PfP'f+PfPj ' 

The momenta k and pj are given by 

fr = x fjJ pl}, p 3 = P» + k^ = 

Note that this construction of momenta is based on the resctriction rrif = 0, which is used 
in the integration of the difference \M. f a ^fjx\ 2 — |-Msub| 2 for mj — > 0. 

In the integration of |.M su b| 2 over the collinear-singular phase space the correct depen- 
dence on a finite mj is required. We sketch this procedure in App. D.2 for a possibly finite 
spectator mass rrij, but here we give only the relevant formulas needed in applications. 
The cross-section contribution a s ^L,fjx that has to be added to the integrated difference 
\M fa ^fx\ 2 ~ \M suh \ 2 is given by 



(5.24) 



(5.25) 



a s f u ^ fX (p f ,p a ;K f ) 



2tt Jo 

where the collinear singularity is contained in the kernels 



dxHf fj>T (P 2 ,x) cr ja ^ jX (k = xp f ,p a ;\ y = TK f ), (5.26) 



?W(P 2 ,a;) = ^ln 



-P 2 (l 



x 



rm 



.f 



P 7f (x)+r(2-x) 



x 



Til 



X) 



X 



(5.27) 
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Figure 5: QED diagrams contributing to e 7 — > e \i fi + at tree level. 

5.4 Phase-space slicing 

Finally, we derive the integral over the collinear phase-space region for the slicing 
approach. This region is defined by restricting the angle O'f between the outgoing and 
incoming / to small values, Q'f < AO <C 1. 

In Section 5.2 this restriction leads to new limits in yfj a , 



,0\2 



<y f ,fa<^(l-x f j a )A0\ (5.2* 



which modify the integrated result to 

H f T f (p° f , x) = \J P °f ( - 1 ~^ Ae ) \ Pf ( x ) + r ( 2 - x )] - - r(l - x), (5.29) 

J V XTTLf / 1 J X 

where the integral is defined analogously to Eq. (5.18). The cross-section contribution for 
the collinear scattering of / is given by 

V C f°adfx(Pf,Va, «/) = / dxHl f (p° f ,x) (T~ /a ^ X (k = Xp f ,p a ; A 7 = TKf). (5.30) 

The same results can be obtained from Section 5.3 with App. D.2, where the new limits 
on the integration variables are given by 



rn 



p 2 



— < *fij < ^2 x m( 1 - *fsj) Ad ■ ( 5 - 31 ) 



6 Application to the process e 7 — > e ^ 

In this section we illustrate the application of the methods described in Sections 3, 
4, and 5 to the process e~7 — > e~/i~/i + at a centre-of-mass energy -y/i much larger than 
the involved particle masses, y/s ^> m e , m^. Of course, this process is not of particular 
importance in particle phenomenology, but it involves the three issues of (i) incoming 
photons splitting into light // pairs, (ii) the collinear production of light // pairs, and 
(iii) forward-scattered fermions and, thus, provides a good test process for these cases. As 
already mentioned in Section 2, our treatment of non-collinear-safe final-state radiation 
has already been tested in other processes. 

To illustrate the formalism, it is sufficient to consider the process e~7 — > e _ /i _ /i + 
in QED, where only the four diagrams shown in Fig. 5 contribute. The corresponding 
helicity amplitudes, including the full dependence on the masses m e and m M , can be 
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obtained from the treatment of e~7 — > e~e~e + presented in Ref. [22] after some obvious 
substitutions. In the following we compare the result with the full mass dependence to 
results obtained with the described subtraction and slicing methods in various kinematical 
situations. Denoting the polar angle of an outgoing particle % by 0j and the angle between 
the two outgoing muons by a^, we distinguish the following cases: 

a) No collinear splittings 

Angular cuts: cut < e - < 180° - cut and M ± < 180° - cut and cut < a^. 
No collinear singularities are included, and the integrated cross section is well defined 
for vanishing fermion masses, i.e. none of the subtraction methods has to be applied. 
The difference between massive and massless calculations indicates the size of the 
fermion mass effects. 

b) Collinear splitting 7 — > e~e + * 

Angular cuts: cut < e - and M ± < 180° — cut and cut < a^. 

The collinear splitting 7 — > e~e + * of the incoming photon is integrated over, so that 

the third diagram of Fig. 5 develops a collinear singularity for backward-scattered 

electrons. The methods of Section 3 are applied to the calculation with massless 

fermions. 

c) Collinear splittings 7 — > / u T /i ± * 

Angular cuts: cut < e - < 180° - cut and cut < a^. 

The collinear splittings 7 — > / u T / u ± * of the incoming photon are integrated over, so 
that the first two diagrams of Fig. 5 develop collinear singularities for backward- 
scattered muons. The methods of Section 3 are applied to the calculation with 
massless fermions. 

d) Collinear splitting 7* — > fi~fi + 

Angular cuts: cut < e - < 180° - cut and M ± < 180° - cut . 

The collinear splitting 7* — > /i~ ji + of an intermediate photon is integrated over, so 
that the last two diagrams of Fig. 5 develop collinear singularities for collinearly 
produced muons. The methods of Section 4 are applied to the calculation with 
massless fermions. 

e) Collinear splitting e~ — > e~7* 

Angular cuts: e - < 180° — cut and M ± < 180° — cut and cut < a^. 

The collinear splitting e~ — > e~7* of the incoming electron is integrated over, so that 

the first two diagrams of Fig. 5 develop collinear singularities for forward-scattered 

electrons. The methods of Section 5 are applied to the calculation with massless 

fermions. 

For the numerical evaluation we set the fermion masses to m e = 0.51099907 MeV and 
m M = 0.10565839 GeV, the fine-structure constant to a = e 2 /(4n) = 1/137.0359895, the 
beam energies to E = E c = = 250 GeV, and the angular cut to 9 cut = 10°. In the 
subtraction and slicing methods the masses m c and m M are neglected everywhere except 
for the mass-singular logarithms, i.e. the laboratory frame defined by the above beam 
energies coincides with the centre-of-mass system. For the fully massive calculation the 
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two frames are connected by a (numerically irrelevant) boost along the beam axis with 
a tiny boost velocity of 0{m 2 / E 2 ). Our numerical results for the different kinematical 
situations and the various methods are collected in Table 1. In addition in Table 2 we 
show the analogous results for the situation where the energy of each final-state lepton 
/ = e - ,/^ is restricted by Ei > 10 GeV. All results are obtained with an integration by 
Vegas [24], using 25 x 10 6 events. While a simple phase-space parametrization is sufficient 
in the subtraction formalism, dedicated phase-space mappings are required to flatten the 
corresponding collinear poles in the slicing approach and when employing the full mass 
dependence of the matrix elements. The fully massive results have been checked with the 
program Whizard [25], where agreement within the integration errors has been found. 

The results obtained with the different subtraction variants, where a spectator is 
chosen from the initial state (IS) or from the final state (FS), are in mutual agreement 
within the integration error, which is indicated in parentheses. Subtraction and slicing 
results are also consistent within the statistical errors as long as the angular slicing cut 
A6 is not chosen too large. For example, some of the slicing results for A9 = 1CT 1 still 
show a significant residual dependence on A9. In the chosen example, the integration 
errors of the subtraction and slicing results are of the same order of magnitude. However, 
we would like to mention that the subtraction approach is often more efficient, as e.g. 
observed in the applications of Refs. [14-17,21] mentioned above. This superiority of the 
subtraction formalism typically deteriorates if complicated phase-space cuts are applied, 
as in the chosen example, because the cuts act differently in the various auxiliary phase 
spaces and thus introduce new peak structures in the integrand. 

Finally, we remark that the impact of mass-suppressed terms is significantly reduced 
if the cut on the lepton energies E\ is applied. This cut guarantees that Ei ^> mi overall 
in phase space, so that mass-suppressed terms are proportional to mf/Q 2 with Q ^ mi. 
Without any restriction on Ei, there are at least small regions of phase space where Q is 
not much smaller than mi, leading to larger mass effects. This feature is clearly visible in 
Tables 1 and 2 when comparing results based on the full mass dependence in the matrix 
elements with the subtraction and slicing results that are based on the asymptotic limit 
mi — > 0. 

7 Summary 

The dipole subtraction formalism for photonic corrections is extended to various 
photon-fermion splittings where the resulting collinear singularities lead to corrections 
that are enhanced by logarithms of small fermion masses m/. Specifically, we have 
considered non-collinear-safe final-state radiation, collinear fermion production from in- 
coming photons, forward-scattered incoming fermions, and collinearly produced fermion- 
antifermion pairs. All formulas needed in applications are provided, only the scattering 
matrix elements for the underlying process and for relevant subprocesses have to be sup- 
plemented in the simple approximation of a massless fermion /. Particle polarization 
is taken care of in all relevant cases, e.g., for incoming fermions and photons. For the 
purpose of cross-checking results in applications, we also provide the formulas needed in 
the phase-space slicing method. 
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Collinear splittings 


Method 


a+_[pbj 


r ii 

a ++ [pbj 


a) none 


full mass dependence 


0.50910(9) 


0.47172(6) 




massless case 


0.51110(9) 


0.47384(7) 


b) 7 — > e~e + * 


full mass dependence 


0.52213(7) 


0.56762(7) 




subtraction (IS spectator) 


0.52424(8) 


0.57027(8) 




subtraction (FS spectator) 


0.52434(7) 


0.57017(9) 




slicing (AO = KT 1 ) 


0.52410(7) 


0.57021(6) 




slicing (A0 = 10~ 3 ) 


0.52431(9) 


0.57021(7) 




slicing (AO = 1(T 5 ) 


0.52423(8) 


0.57028(7) 


c) 7 -> ^ T /^ ± * 


full mass dependence 


2.5890(5) 


2.3615(4) 




subtraction (IS spectator) 


2.5872(3) 


2.3586(5) 




subtraction (FS spectator) 


2.5873(8) 


2.3585(5) 




slicing (A0 = 10" 1 ) 


2.5883(3) 


2.3609(2) 




slicing (AO = 10~ 3 ) 


2.5859(8) 


2.3578(8) 




slicing (A6 1 = 10 5 ) 


2.5876(13) 


2.3572(13) 


d) 7* -> 


full mass dependence 


0.54076(8) 


0.53357(8) 




subtraction (IS spectator) 


0.54309(8) 


0.53597(7) 




subtraction (FS spectator) 


0.54306(8) 


0.53603(7) 




slicing (AO = KT 1 ) 


0.53164(19) 


0.52386(16) 




slicing (AO = KT 3 ) 


0.54287(17) 


0.53624(15) 




slicing (AO = KT 5 ) 


0.54335(18) 


0.53580(18) 


e) e~ — > e"7* 


full mass dependence 


5.5465(7) 


4.7060(6) 




subtraction (IS spectator) 


5.5495(4) 


4.7070(3) 




subtraction (FS spectator) 


5.5484(6) 


4.7064(5) 




slicing (AO = KT 1 ) 


5.5313(1) 


4.6880(1) 




slicing (AO = 10" 3 ) 


5.5488(3) 


4.7071(3) 




slicing (AO = KT 5 ) 


5.5486(5) 


4.7067(4) 



Table 1: QED cross sections o n \ for e~7 — > e~/i~/i + in the various setups described in the 
main text, with signs k and A of the helicities of the incoming e~ and 7, respectively. The 
results are obtained with the indicated methods, where IS and FS stand for spectators in 
the initial and final states, respectively. 
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Collinear splittings Method (j + _[pb] a ++ [pb] 
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9 zLQ9^(m^ 






slirinp- (AO — 1 0~M 

OllLjlllti I Z_lL/ 1U J 


9 zlnQ^f^ 

Zv . t:(J UxJ I Zv J 


2 2606f2"l 




slirina- ( AO — 1 fl -3 ^ 

OllLlllg, \I—\U — 1U 1 




9 9fi98('7'l 




slirinp (AO — 1 O -5 ^ 

OllL/lllIZ 1 i—WJ 1U ) 


2 490^(1 9"l 


9 9696f 1 2 1 ) 


clj 7 


full mass dependence 


U.4o0U0^ ( ) 


U.4 i tSau^o J 




oUU 11 dL L1U11 I lu opLLldLiUl 1 


fl 4R<S9PY7 N I 






subtraction (FS spectator) 


U.4o0oU^0 ) 


U.4 I 4U1 ^0 J 




slicing (AO = 10" 1 ) 


0.47588(19) 


0.46363(13) 




slicing (A6> = 1(T 3 ) 


0.48607(19) 


0.47399(14) 




slicing (AO = 1(T 5 ) 


0.48623(20) 


0.47425(15) 


e) e~ — > e"7* 


full mass dependence 


5.4878(6) 


4.6467(5) 




subtraction (IS spectator) 


5.4866(3) 


4.6471(3) 




subtraction (FS spectator) 


5.4871(5) 


4.6475(5) 




slicing (AO = KT 1 ) 


5.4690(1) 


4.6278(1) 




slicing (AO = 1(T 3 ) 


5.4869(3) 


4.6467(3) 




slicing (AO = 1(T 5 ) 


5.4862(5) 


4.6466(4) 



Table 2: Same as in Table 1, but with an energy cut E\ > 10 GeV for all final-state 
leptons. 
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As an example illustrating the use and performance of the proposed methods we have 
explicitly applied the subtraction procedures to the process e~7 — > e~fi~fi + and compared 
the results to those obtained with phase-space slicing. The presented subtraction variants 
will certainly be used in several precision calculations needed for present and future collider 
experiments such as the LHC or ILC. 
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Appendix 

A More details on non-collinear-safe final-state radiation 



Here we generalize the results of Section 2.1, where non-collinear-safe photon radiation 
off fermions is treated, to the situation where massive spectators in the final state exist. To 
this end, we only have to consider the case of final-state emitter and final-state spectator. 

For rrii — > 0, m 7 = 0, but rrij ^ 0, the boundary of the yij integration [given for the 
massless case in Eq. (2.8)] is given by 



yi(z) = 



m 2 (l-z) 



PfjZ 



i{z) + 1 + ^(z)Uz) + 2] 



with £(z) 



mi 



2P 2 j z{\ 



(A.l) 



and the functions relevant for the integrand gfj^f 1 behave as 



Pik 



Rij(y) = 



\ 



(i-vY 



4mjy 



P 2 ' 



l - 



2m|^ 

p?. i- 



(A.2) 



The evaluation of Eq. (2.6) now becomes non-trivial and yields 



e$$\p%,z) = -Pff(z)Hi£2[i-v{z 



zP 2 - 
13 



+ (l + z) ln(l + 



m\ 



+ (1 + z) ln(l - z) - 

2 Wi+ 



V PMz)) (l-z)a(z) 



2z 
1-z 

P 2 V (z)[l - z V (z)\ 
m 2 Al — z) 



-21nfl-i^ 

V 1 + a(z) 



+ a(z) In 



'l-z))}-Q^{^ v zl 



gj£\p*,z) = i-z, 



(A.3) 
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with the auxiliary functions 

a(z) 



4m f / 1 \ i \ I [1-2/2(2)12 for 2; < I, . A 

X 1 + ^/2(1-2), V{z) = \ l h ' w ? (A.4) 

\J I [l-y 2 (*)](l-,z) for z > - 2 . 



For mj — > 0, the results for ^^(P 2 -, 2^) reduce to Eq. (2.10), as can be easily seen after 
realizing that r)(z) = 0(m,j) and a(z) = 1 + 0(m?) in this limit. 

B More details on the subtraction for 7 — > ff* splittings 

B.l Factorization in the collinear limit 

In this section we derive the asymptotic behaviour (3.3) of the squared amplitude 
\M. ia -*fx\ 2 for the case where the outgoing light fermion flies along the direction of the 
incoming photon. We consider polarized incoming photons with momentum /c M and po- 
larization vector , where A 7 = ± is the sign of its helicity. We further introduce a 
light-like gauge vector (n 2 = 0,nk 7^ 0), i.e. is characterized by 

£^ = ( £ y*, ke Xl = ne Xl = 0. (B.l) 

In the following we make use of the identity 6 

where 



4(4)* = 44 = ^(k) T ^ VP °K^ (B.2) 



hV^v _|_ rjMi^ 

E^{k) = e\e v _ + ete\ = -<T + — — — — (B.3) 

kn 

is the polarization sum of the photon in four space-time dimensions and e^ vpa the Levi- 
Civita tensor with e 0123 = +1. 

In a gauge for the photon where nk = 0(k°), it is easily shown by power counting that 
the logarithmic singularity arising from the phase-space region kpf — (D(mj) (m/ <C k°) 
originates from diagrams in which the incoming photon collinearly splits into a light //* 
pair. The generic form of such graphs is shown in Fig. 6. Assuming summation over the 
polarization of the outgoing fermion /, the squared matrix element, thus, behaves like 

\Mja-*fx(k,p a ,Pf, A 7 )| 2 

J^O QV f fa^X (Pf, Pa) ~j / + 7 ^ 7 (h + m f )^ (Pf, Pa) , (B.4) 

1 Pj — / Pf — / 



6 This identity is easily proven using a representation of the polarization vectors by Weyl spinors. 
Employing the conventions of Ref. [23], we have e AB = £+& AB = V2n A k B /(kn) and e AB = e^L<T AB = 
\/2k A n B I (kn)* for the polarization bispinors, so that 
e^k p n a = i(e A6 e dd e BD e FH - e A ^e^e BF e DH )^ AB ^ dD a^ dH k p n a 
= \(k A n^e BD k x n x - e Ad k B n D k x n x )a\ B a v dD 

= H_ k A n B n C k D + n A k B k C n D^^^ 

= ±(kn)(kn)*{e> 1 + e v _ - ete v + ) = i(kn)(s^s v _ - ete v + ). 
The only non-trivial step is the third equality which follows from a twofold application of Schouten's 
identity. 
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a 1 a 1 

Figure 6: Generic squared diagram for the splitting 7 — > f f* (left) and the corresponding 
squared diagram for the related process with an incoming / (right). 



where Tj a ^ x (pj,p a ) includes all information of the subamplitude indicated by the open 
blob in Fig. 6 and Tj a ^ x = {Tj a ^ x )^Q. To leading order in m/ — > 0, the squared 
amplitude for the subprocess fa—>X can be written as 

\Mf a ^ x (pf, p a ; Kf)\ 2 = T f - a ^ x {p f -, p a ) [w- Kf tf + 0(m f )] T fa ^ x (pf, p a ), (B.5) 

with u± — |(1 ± 75) and kj = ± denoting the sign of the / helicity. In order to find the 
relation between these squared matrix elements, we insert identity (B.2) into Eq. (B.4) 
and eliminate the e tensor via Chisholm's identity, 

■ l( nW ls = (7 * 7 /y _ g ^ + g *nj> _ /-y*) T5) (B .6) 

i.e. we trade the e contributions for a 75 insertion in the spinor chain. Next, we isolate 
the leading terms in the collinear limit kpf = 0{nv*) — > 0. This limit can, e.g., be 
parametrized by the decomposition of the momentum of / 

p» = (l- x )W+p^+pl r (B.7) 

with x = 1 — p°f/k°, kpf^x = 0, and pf )V = (where boldface symbols refer the spatial 
parts of momenta). In this decomposition we have O(p f ± ) = O{p f r ) = 0(mj) and 
P/_i_ = 0(mj). Thus, each component of the orthogonal 3- vector p^j_ is of 0(rrif). After 
some straightforward simplifying algebra, the result of applying the power counting to 
\M ia ^ fX \ 2 is 7 

\M 7a ->fx(k,Pa,Pf] A 7 )| 2 



kp f ^o 2x(kpf) 



Tfa^ X (Pf,Pa) \ 1-2X(1-X) + -I- 1 j)J 



kp f 



7 Actually there are also terms proportional to A 7 m// (kpf) 2 2y o _ >x -^/ ) x7sTj _ > ^- , which at first sight 
seem to contribute in 0(mj 2 ) in the limit mj — * 0. Although these terms obviously disappear from the 
subtraction function after setting rat to zero, they potentially contribute to the corresponding integrated 
subtraction terms, in which the limit mf — > is taken after the singular phase-space integration. However, 
the integration over the azimuthal angle of p/, which is always assumed in our analysis, leads to a further 
suppression by one power of m/, so that the contribution to the phase-space integral of \A4 ia ^fx\ 2 
is mass suppressed. Thus, these terms are irrelevant in the construction of a subtraction function to 
separate mass-singular terms in the collinear cone. 
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-A, 



2x-l + 



m 



kp f 



l4f\ T ja^x{Pj,Pa) 



1 



kl ^o 2 Q " e2 \( h + + ) [\ M f^x(Pf,P a ; +)| 2 + |A4/«^(p/,p«; -) 

+ Xj(h1 f - hl f ) \\Mf„ x (pj,Pa t +)| 2 - \M fa ^ x (p f -,p a ; -)f 



(B.8) 



The last form results from the last but one by simply substituting \Mf a _, x \ 2 and the 
K± functions defined in Eq. (3.4), whose arguments are suppressed in the notation. This 
completes our proof of Eq. (3.3), which is a more compact version of this result. 

In this section we have explicitly treated / as fermion and / as antifermion. The 
opposite case with / being an antifermion and / a fermion is obtained analogously and 
leads to the identical final result (3.3), although some signs in intermediate results are 
different. This fact is, of course, to be expected, because relations between squared helicity 
amplitudes cannot depend on our convention which fermion we call the antiparticle of the 
other. 

B.2 Dipole subtraction for 7 — > //* splittings with massive final-state spec- 
tator 

Here we give some details on the derivation of the integrated subtraction part presented 
in Section 3.3 for the collinear splitting 7 — > //* in the process 7a — > fjX, where j is a 
possibly massive spectator. We start by generalizing the form (3.24) of the new momenta 
upon restoring the correct dependence on m/, 



P 2 + m 2 -m 



/v-v _ p2 y° 1 2 p2- j 2 P 2 
Pj ( x ) = + Pf(x) , p] = p 1 * (x fj „), 
where the following shorthands are used, 



P 2 -m 2 



> A/ il7 = A(P 2 ,mJ,m 2 ), 



with the auxiliary function 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 



(B.9) 
(B.10) 
(B.H) 



The new momenta satisfy the on-shell conditions p 2 j = m 2 , pj = m 2 and correctly behave 
in the collinear limit, pj — > xk, where kpf = 0(m 2 ) — > 0. The splitting of the (N + 1)- 
particle phase space into the corresponding iV-particle phase space and the integral over 
the remaining singular degrees of freedom is given by 



d<t>(Pf,Pj,k x ;k+p a ) = J dx j d(j){p j (x),kx;Pf(x)+p a ) J[dp f (P 2 ,x, z fjjl )\, (B.12) 
with the explicit parametrization 
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and kx denoting the outgoing total momentum of X. The upper kinematical limit of the 
parameter x = Xfj a is given by 

_p2 

Xl = — =" ; . (B.14) 

-P 2 + 2m f m j K J 

The integration of the azimuthal angle 0/ of / simply yields a factor 2n, but the integration 
of the auxiliary parameter 

z . _ kpj ,g ^ 

fhy kpf + kpj 1 ' ; 

with the boundary 

2m 2 x + P 2 (x - 1) =F \/P 4 (l - x) 2 - Am 2 f m 2 x 2 

- - — 2 U-U — — < B16 > 

is non-trivial. The integration kernels occurring in the final result (3.25) are defined as 

and can be evaluated without problems analytically (even for finite to/) yielding Eq. (3.26) 
for rrif — > 0. 



= r™dz fj „hlf T (k,p f , Pj ) (B.17) 



C More details on the subtraction for 7* — > // splittings 

In this appendix we supplement Section 4.3, where the subtraction procedure for 
collinear 7* — > // splittings has been described for a final-state spectator j. In the 
following we fully take into account the spectator mass m,j. The derivation widely follows 
Ref . [8] , where the treatment of the g* — > QQ splitting with a massive quark Q has been 
considered. Our approach differs from the one described in Ref. [8] in the level of inclu- 
siveness that is assumed in the collinear limit; in contrast to that paper we do not assume 
a recombination of the // pair in the collinear limit, but instead control the individual 
momentum flow of / and /. 

For arbitrary mass values to/ and rrij the subtraction function can be constructed as 
in Eq. (4.17) with the generalized radiator function 

,2 



^(p/,P/,Pi) = (p/+p/)2 „. 



(Pf + Pf) 



\ Z ffj P f Z ffjPf) \ Z ffj P f Z ffj P f) 



■(CI) 



In addition to the parameters t/ffj and Zjfj, which are defined as in Eq. (4.19), we make 
use of the following auxiliary quantities, 

P 2 = p 2 - 2m) - , 

P 2 Vm - 2mJ _ ^[2m 2 + P 2 (l-y ffj )] 2 -4m 2 P 2 



v f 



\P 2 y ff - j + 2my V * P^-Vffl) 
-(IT VjVf), z) f > = Zffj - -(1 - Vj), z) f > 



zi, 2 — ^(1 T VjVf), z f f = Zffj - -(1 - Uj), 4?] = *ffi - ^ - ( C - 2 ) 
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The parameter k is arbitrary, because the singular behaviour does not depend on it; in 
practice the independence of the final result on k can be used as check. The auxiliary 
momenta entering the hard scattering matrix element for the subprocess ab — » 7j'X also 
become more complicated, 



P 2 - m) 



P 



M _ \ _|_ 



^\(P 2 ,(p f +p f ) 2 ,m 2 )V J P 2 



P2 + m i p. 
2P 2 



(C.3) 



In order to integrate the subtraction function we need the azimuthal-averaged version of 

h 7f,P 



h f fj(pf,Pf,Pj) 



(Pf+PfYvj 



Pf 1 {z f!] ) + 2{l-K)z 1 z 2 + 



2nm 2 



(C.4) 



(Pf+Pf) 2 

and an appropriate splitting of the phase space of the momenta pf, pj, pj, 

J d4>(p f ,p f , Pj ;P) = J dtfapfiP) J [dp f (P 2 ,y fIj ,z ffj )], (C.5) 

J [dPf{p ' y ^ z ^ )] = M^y L dVffj (1 " Vffl) w,) dZffl J d0/ ' 

where 



Vi = 



2mJ 

p2 ' 



2/2 



2m j (VP 2 -m j ) 
P 2 



(C.6) 



and z lt2 (yffj) are the Z\^ of Eq. (C.2), evaluated as functions of Vffj- Up to this point, 
the full dependence on mj and rrij is kept. 

Since we want to keep the momentum flow in the collinear limit open, i.e. the Zffj 
integration should be done numerically, we have to interchange the order of y^jj and Zfjj 
integrations in the singular phase-space integration over /[dpj]. For arbitrary masses rrif 
and rrij, this seems hardly possible analytically, so that we focus on the limit mj — > in 
the following, because this is the interesting case. We define 



p* rvi{z) 

n ffA P i z ) = -7T / (S dVffjV-VffjjhffjiPfiPfiPjh 

£ J'm(z) 



2 Jyi(z) 



HffA p2 )= f'dzHfjj^z), 



(C.7) 



where we were allowed to use rrif = in the prefactors and in the integration limits of 
z = Zjjj. The relevant asymptotics of yi^( z ) for rrif — > is 



m 2 f z 2 + (1 - zf 



V2{z) 



AP 2 z(l 



rrij — rrij 



P 2 z(l - z) ' 
The actual integration over y^jj yields 

' ^4P 2 z(l- z) +m 2 - 



^jAP 2 z(l - z)+m 2 + rrij 



(C.l 



n ff j(P i ,z) = p fJ (z) 



2 In 



2m f 



~ 1 - V(z) 
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-21n[l-77(*)] + 



m 2 [l — rj(z)] 
m 2 + rj(z)P 2 



2ui 2 / P 2 \ 

+ -^(l - K + z 2 + (1 - ^) 2 ) ln(l + 77(^-2 ) + 2z(l - z) (C.9) 



with 



[ [1 -y 2 (^)](l - z) for ^ > i. 

The case rrij = given in Eq. (4.22) can be easily read off after realizing that rj(z) = 
0(rrij). For the evaluation of Hffj(P 2 ) it is easier to integrate first over z and then over 
Uffj- The result is 

tt 4 , fVP 2 -™^ 16 Am j ( 2\2m 2 ( 2m j 
HffAP ) = - hi 3 - 1 =1 h (« Hf n ; 

(C.11) 

which could also be derived from Eq. (5.36) of Ref. [8]. For rrij = this obviously leads 
to the form given in Eq. (4.22). 

D More details on the subtraction for / — > fj* splittings 

D.l Factorization in the collinear limit 

In this section we derive the asymptotic behaviour (5.2) of the squared amplitude 
\Aif a ^fx\ 2 for the case where the incoming and outgoing light fermions become collinear. 
We consider polarized incoming fermions / with momentum p^ and helicity of sign Kf — ±. 
The corresponding Dirac spinor u(j)f, Kf) is an eigenspinor of the helicity projector 

S «/ = ^(1 + k/7s/* P/ ), (D.l) 

where the polarization vector 

3 » (MAe,) (D. 2) 

is aligned to the direction e/ = p//|p/| for helicity eigenstates. Defining the light-like 
vectors /c M = fco(l,e/) and = (1, — e/), the polarization vector s£ can be decomposed 

into k^ and as follows, 

rf, = p^~k» - -pL^n». (D.3) 
p f 2m f k 2{ Pf n) v ; 

Note that the momentum k^ of the virtual photon fulfills kn = O(k ) in the collinear 
limit, because then /c M = /c M + Oimf). The vector n M will be used as gauge vector in the 
explicit definition of photon polarization vectors for the subprocess 7a — > X below. 

Power counting reveals that the logarithmic singularity arising from the phase-space 
region Pfp'f = 0(m 2 ) — > (rrif <pj) originates from the square of diagrams in which 
the incoming fermion collinearly emits a photon that triggers the production of X. The 
generic form of such graphs is shown in Fig. 7. Assuming summation over the polarization 
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a 1 a 1 

Figure 7: Generic squared diagram for the splitting / — > fj* (left) and the corresponding 
squared diagram for the related process with an incoming 7 (right). 



of the outgoing fermion /, the squared matrix element behaves like 

\Mf a ^fx(Pf,Pa,P'f;K f )\ 2 

N 2 e 2 

pfir° ^{^f + m fhM + m f H} T^ x (k,p a )*T^ x (k,p a ), (D.4) 

where T^ a ^ x (k,p a ) includes all information of the subamplitude indicated by the open 
blob in Fig. 7. To leading order in mj- -> 0, the squared amplitude for the subprocess 
7a — > X can be written as 

\M ja ^ x Ck,Pa,\,)\ 2 = el^Ka^xC^PaYex^T^ikpa), (D.5) 

with the helicity A 7 = ± of the incoming photon and the light-like vector In order to 
relate the fa process with the 7a subprocess, we first evaluate the trace in Eq. (D.4) and 
drop all terms that vanish owing to the Ward identity k^T^ a ^ x (k,p a ) = 0. Inserting the 
form (D.3) of s£ , the result can be written as 

\Mf a ^fx(Pf,Pa,p'f]K f )\ 2 

_ N c , f Q}e 2 f A PftllPftV m f ba f(Pf n )l0 m ) r 0\ 

PfP'f^ -k* I 9ilV k 2 k 2 6 ^0 K { k() fan) J. 

xT^ x (k,p a )*T^ x (k,p a ). (D.6) 

Now we make use of the collinear limit which is characterized by pjp'f = m 2 — pfk = 
0(m 2 ) — > 0. We decompose the photon momentum according to 

fc" = X p) + k^_ + k? (D.7) 

with x = k°/p°f, Pfk± = 0, and k r = 0. In this decomposition we have O{k L ) = 0(k®) = 
0(m 2 ) and = (D(m 2 ), i.e. the vector kj_ can be counted as 0(rrif). Moreover, we get 
k\ = x 2 m 2 + k 2 (l — x) + 0(m^). In the determination of the leading collinear behaviour of 

Eq. (D.6), we can replace the momentum fc M by the light-like momentum k^ = k^ + Ofaf) 
in the two T ja ^ x (k,p a ) terms. The expansion of the two terms with the e-tensor is also 
straightforward. With the help of identity (B.2), the contraction t ilva pk a, nP becomes 

Va^ a ^ = W/3^" n/3 + 0(m f ) = i(kn) e + ^(k)e^ M (k) - £^ tfl (k)e +tl/ (k) + 0(m f ). 

(D.8) 
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The second contraction z^apk^kP can be expanded upon writing e^ vol p = 
9^ 9v v 9a a Cu'v'a'p with the following decomposition of the metric tensor, 

g ^ = -L^V + Pra") - - e^(k)e v + (k). (D.9) 

2/c 

The e-tensor now only appears as e flua ^e'^(k)e l L(k)n a k 13 = 2i/c , and the momentum k^ 
with an open index can be replaced via 



2ko_ 
(kn) 



fc" + (e_(jfe) • fc) 4(jfe) + (e+(jfe) • k) e^(jfe)] + O(mJ), (D.10) 



which follows from Eq. (D.9) upon contraction with k u . This procedure spans the tensor 
e liua pk a k 13 in terms of e± t ^(k)e TtU (k) and covariants involving /c M or k v . The latter do not 
contribute because of the Ward identity /c M T^ a ^ x (A;,p a ) = 0. The expansion of the various 
scalar products for mj — > is straightforward, yielding 

~ B ik [k 2 (x — 2) — x 2 m 2 f ] r ~ ~ ~ ~ -, 

e^ a /3& fc p = xjp~n) *A&) - 

+ (terms proportional to fc M or k u ) + (9(my). (D.ll) 

Inserting Eqs. (D.8) and (D.ll) into Eq. (D.6) and performing the expansion in the 
collinear limit leads to the form given in Eqs. (5.2) and (5.3). 

The final step of performing the azimuthal average around the collinear axis, which 
leads to Eqs. (5.5) and (5.6), is most easily carried out by fixing a specific coordinate 
frame. In a frame where the direction of p/ is given by ej = (0, 0, 1), the vectors k^ and 
e±(k) are given by 

k" = (fc , 0, 0, fc ), 4(*) = ^(0, 1, ±i, 0). (D.12) 

Recall that and differ only by terms of 0(m,f) in the collinear limit. According to 
definition (D.7) the leading term of takes the form 

k± = (O,-|k ± |cos0 / / ,-|k ± |sin0 / / ,O) + O(mJ), (D.13) 

where 0j is the azimuthal angle of p'j = p/ — k. In this parametrization the average 
(kj_k^) t// is easily calculated to 

<*MV, = - Ydiag(0, 1, 1, 0) + 0{m)) = -^E^(~k) + 0(m% (D.14) 

while it is trivially seen that the tensors e±(k)e±(k)* do not change after taking the 
azimuthal average. With these considerations the transition from Eqs. (5.2) and (5.3) to 
the averaged form in Eqs. (5.5) and (5.6) is straightforward. 
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D.2 Dipole subtraction for / — > fj* splittings with massive final-state spec- 
tator 

In Section 5.3 we have presented all formulas needed for the case of a massless final- 
state spectator in practice, but did not go into the details of their derivation. Here 
we close this gap by deriving the formalism in the more general situation of a possibly 
massive spectator j. We keep the general definition (5.22) of the subtraction function, 
but generalize the subtraction kernel as follows, 



4(1 - x) k^k u L (p f - p' f ) 2 (l -x) + m)x 2 



-gT- 



+ 



11 

X 



2-x + 



2x 2 m 2 
(Pf ~ P'f) 



k\ (p f - P ' f ni-x) 

e^{k)*e v + (k) - e»(k)*e u _(k)) 

(D.15) 



because we need the correct dependence on the emitter mass to/ for the integration of 
|-M S ub| 2 below. The auxiliary parameters still have the form (5.24), but the auxiliary 
momenta become more complicated, 



m 2 ~P 2 



x 



-P 2 R(x) 
Pj(x) = P M + F(x), 



pj + 



P 2 + 2m 2 f x 



2xP 2 
Pi = Pj( x fi,f) 



TO' 



P 2 



2P 2 



-P", fr* = Hx fj , f ), 



Pjk ,,, p'jk 

k± = ttP/ ^jPy 
Pjk pjk 

Here we made use of the abbreviations 



P^ = p'f + 



P 2 = P 2 - 2m) 



TO 2 , 



R(x) = 



(D.16) 



^(P 2 + 2to 2 x) 2 - 4x 2 m 2 P 2 
Zpl 



(D-17) 

The new momenta satisfy the on-shell conditions k 2 = 0, p 2 = m 2 and correctly behave in 

the collinear limit, k — > xpf, where p/p'j = 0(m 2 ) — > 0. Before carrying out the singular 
integrations, we average the subtraction function over 0j, yielding 

(|M,ub(«/)lV = N c , f Q)e 2 h^( Pf ,p' f ,p 3 ) \M ia ^ jX (k,p a ; A 7 = tk s )\ 2 (D.18) 



with summation over r = ± and 



R ff( x fjj) + 



2m } x kf 



/,2 "-,o./ + -"'/•'•//./ 

+ r\2- x fjJ + 



2x %,f m } 



/rJ "-//,/ • -'"Vfj./J. 



(D.19) 



The splitting of the (N + l)-particle phase space into the corresponding iV-particle 
phase space and the integral over the remaining singular degrees of freedom is given by 

J d<P{p'f,Pj,kx;Pf+Pa) = J dx J d(f)(p j (x),k x ;k(x)+p a ^ J[dp' f (P 2 ,x,z fjJ )}, (D.20) 
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with the explicit parametrization 



1 s —P 2 rM x ) 
A(2~nfl x 2 R(x) J Z1 ( X ) dZfj ' J 



(D.21) 



The upper kinematical limit of the parameter x = Xfjj is given by 

-P 2 

Xi 



—P 2 + Imjvrij 



(D.22) 



The integration of the azimuthal angle 0j of f(p'f) simply yields a factor 2ir. The non- 
trivial integration over Zfjj has the boundary 



z li2 {x) 



2mjx + P 2 (x - 1) =f R(x)JP 4 (l - x) 2 - Am)m 2 x 2 



2[P 2 (x- 1) + (m? f + m 2 )x] 
Defining the integrated subtraction kernel according to 

_p2 . Z2 ( x ) 



(D.23) 



n f j f T {P 2 , X) = dZfjJ h f jf T (p f , P ' f , P j), 



(D.24) 



the cross-section contribution of the subtraction part takes the form (5.26) in the limit 
rrif — > 0. For a non-zero spectator mass m,j, the function TCjf T (P 2 ,x) reads 



In 



P 4 (l 



2 V x 3 m} [-pa ( i _ x ) + m 2 x j 



P lf (x)+r(2-x) 



1-x 



x 



-t(1-x). 
(D.25) 
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